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Abstract. Molecular transitions recently discovered at redshift z-^bs = 2.059 toward the 
bright background quasar J2123-0050 are analysed to limit cosmological variation in 
the proton-to-electron mass ratio, fi = m p /m s . Observed with the Keck telescope, the 
optical spectrum has the highest resolving power and largest number (86) of H 2 tran- 
sitions in such analyses so far. Also, (7) HD transitions are used for the first time to 
constrain yu-variation. These factors, and an analysis employing the fewest possible free 
parameters, strongly constrain yu's relative deviation from the current laboratory value: 
Ayu/yU = (+5.6 ± 5.5 sta t ± 2.7 sys ) x 10~ 6 . This is the first Keck result to complement recent 
constraints from three systems at z abs > 2.5 observed with the Very Large Telescope. 
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1. Introduction 

The assumed invariance of fundamental con- 
stants within the Standard Model of parti- 
cle physics must be tested experimentally. 
Temporal or spatial deviations from the cur- 
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rently measured laboratory values of impor- 
tant constants, such as the fine structure con- 
stant, a, and the proton-to-electron mass ratio, 
/j = m p /m e , could point to more fundamental 
theories of physical interactions. 

Limits on the variation of fi have been 
obtained in Earth-bound lab oratory exper- 
iments (e.g. iBlatt et all I2008I) but can also 
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be derived over larger cosmological scales 
from molecular hydrog en absorpt i on in 
the spectra of quasars (fThompsonl [1971 



Varsh alovich &LevshakovL Il993l) . Variation 
of ju manifests in a mass-dependent velocity 
shift, Av,-, of a ro-vibronic transition i. The 
magnitude and direction of Av, is characterised 
by a sensitivity coefficient Kj. Measurement of 
a shift in the observed redshift of the transition, 
Azi, relative to the redshift of the absorption 
cloud, z a bs, in which the transition is observed, 
implies a variation of yiz, 



Av,- AZi A/J. 

— ~ &i — , 

C 1 + Zabs V 



(1) 



where A/u/fi = (fi z -/WMab fbr//i ab , the cur- 
rent laboratory value of fi and fi z the value in 
the absorption cloud. Using two or more tran- 
sitions allows for Aju/// and z a bs to be simulta- 
neously measured. 

Current studies derive limits on Ayu/// us- 
ing a small sample of H2 absorption systems 
from the UVES instrument on the Very Large 
Telescope (VLT) in Chi le (llvanchik et all 
l2005tlReinhold et all l2006). More recent anal- 
yses, also VLT-ba s ed, yield null con s traints 
king et all l2008t IWendt & ReimersL l2008t 



iThompson et all 120091) . 

In this paper we describe our analysis of a 
new H2 & HD absorber observed with the Keck 
telescope. F ull details of the a nalysis will be 
presented in lMalec etaTI (l2010h . 
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Fig. 1. Some of the 86 H 2 and 7 HD lines from the 
J2123-0050 Keck spectrum on a velocity scale cen- 
tred at z abs = 2.0594. The molecular lines fall in the 
wavelength range 3100-3421 A with SNR 7-25 per 
1.3kms~' pixel. The spectrum (black histogram) 
is normalized by a nominal continuum (upper dot- 
ted line). Also shown are the fiducial 4 component 
fit (solid curve, VCs marked by tick marks), local 
continua and zero levels (upper and lower dashed 
lines) and hya lines (denoted by tick marks offset 
above the molecular VCs). Note that the two left- 
most components are nearly coincident in velocity. 



2. Data 

2.1. Keck spectrum of J21 23-0050 

We study the newly discovered H2 absorber 
at Zabs = 2.059 towards the z em = 2.261 
quasar SDSS J212329.46-005052.9 (hereafter 
J2 123 -0050; Milutinovic et al., in prepara- 
tion). Given the scarci ty of such systems 
(iNoter daeme et al., 2008) and further selection 
criteria imposed by varying-yU analyses (z^s > 
2 to shift enough transitions above the atmo- 
spheric cutoff, high enough H2 column densi- 
ties and background quasar luminosities) this 
new absorber is exceptional. The R m 110000 
spectrum was obtained using the HIRES in- 
strument on the Keck telescope in Hawaii and 



covers wavelengths 307 1-5896 A. We find 86 
H2 and, for the first time in such analyses, 7 
HD transitions usable for constraining Afi/fi. 
Notably, this is also the first time Keck/HIRES- 
based data is used to constrain //-variation with 
similar precision to the previous VLT results. 

Figure Q] shows a small, but representa- 
tive subset of H2 and HD transitions in the 
data, with a range of overall line-strengths, 
SNRs and ground-state rotational levels, char- 
acterised by the quantum number J ('/- 
levels'). Two distinct spectral features (SFs) 
are seen, separated by ^20kms~'. This is dis- 
tinct from the velocity components (VCs) used 
to model the molecular absorption system. In 
Section 13.21 we find evidence for 4 VCs given 



884 



Malec: Limit on the variation of the proton-to-electron mass ratio 



the SNR of the data. The strong left-hand SF 
appears saturated for most low-/ transitions 
and the weaker right-hand SF appears unsat- 
urated in almost all transitions. Only the left- 
hand SF is detected in the HD transitions. 



2.2. Molecular data 

Quasar-based varying-yU analyses rely on sub- 
pixel measurements of deviations of astro- 
nomical transition wavelengths from labora- 
tory wavelengths, and on the sensitivity coeffi- 
cients, Kj, that characterise the deviations. For 
H2, we use the most accu rate laboratory wave - 
lengths avail a ble fro m lAbgrall et al. ( 1993h. 
Ubac hs et alJ ([2007) and ISalumbide s et alJ 
(2008). For HP , we u se the wavelengths listed 
in llvanovetaT] (120081) . The H 2 a nd HD Kj co- 
efficie nts u sed were calculated by lUbachs et alJ 
d2007l) and llvanov et all d2008l) . respectively. 



3. Analysis and Results 

3.1 . Fit and x 2 minimization analysis 



Previous analyses (e.g . Ilvanchik et all l2Q05b 
iThompson etail [2009) 

used a 'line-by line' 
fitting approach. The redshifts of individual H2 
lines were derived from independent fits. Using 
a linear fit to the reduced redshifts versus the K, 
coefficients, the value of A/i//i can be derived 
as defined in equation This approach re- 
quires that only a single VC is fit to the molecu- 
lar lines. Lines blended with neighbouring fea- 
tures were excluded from these analyses. 

All of the molecular lines of interest in our 
data fall in the Lya forest. Many are blended 
(to varying degrees) with Lyo- absorption lines. 
Some are blended with metal lines in the same 
absorption system (we model this absorption 
accurately using transitions that fall outside of 
the Lya forest). The molecular absorption it- 
self shows two SFs which, while distinct, are 
'blended' . To utilise the many H2 lines present 
in the data, a number of which would other- 
wise be rejected, H2 and HD absorption pro- 
files, together with the broader Lyman-o- pro- 
files, are fit simultaneously. With a lesser influ- 
ence on the derived value of Afi/fi, there also 
exist local uncertainties in the zero and con- 



tinuum flux levels, set as free parameters in 
the fit. This 'simultaneous fitting' technique al- 
lows Ajtz/yU to be calculated with the maximum 
precision available in the data, while taking 
into account the uncertainties brought about 
by ble nds, level adjustm ents and also varying 
SNRs dKing et all I2OO8L describes further ad- 
vantages of this technique). 

Each line, or VC, associated with an ab- 
sorbing cloud, including Lyo- and metal lines, 
is modelled as a Voigt profile, parametrized by 
the Doppler width, b, column density, N, and 
redshift, tabs- Different molecular transitions 
share these parameters in physically meaning- 
ful ways. For H2 and HD we assume a common 
velocity structure: a given VC shares the same 
Zabs m a ll transitions. Each VC has a common 
N in the same /-level and a common b for all 
/-levels. Because HD lines are few and rela- 
tively weak the HD N values are constrained 
to follow the same pattern across the profile 
as the H2 /=0 values. In relating the absorp- 
tion cloud parameters the number of free pa- 
rameters in the absorption model is minimized. 
We relax these assumptions as part of the inter- 
nal consistency checks described in Section[33| 
and find no evidence of strong deviations from 
our final result. The parameter of key interest, 
Aju/yU, is common to all molecular VCs. 

The x 1 between the data and the absorp- 
tion model is minimised using the program vp- 
fit, designed for modelling complex, interre- 
lated absorption lines with the ability to link 
physically related parameters as defined in our 
model. Statistical 1-cr uncertainties on best fit- 
ting parameters, including Ayu/yU, are derived 
by vpfit from the final parameter covariance 
matrix. The reliability of the calculated uncer- 
tainty in A/i/yU is tested using fits to simulated 
spectra in Section POl 



3.2. Fiducial absorption model 

The initial model for the molecular absorption 
used two VCs, one for each SF. The Lya for- 
est lines, blending metal lines, continuum and 
zero levels were added and refined assuming 
this preliminary molecular model. The veloc- 
ity structure was progressively built by adding 
molecular VCs. After each fit was optimised, 
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Fig. 2. Composite residual spectra (CRS), formed from 23 relatively unblended H 2 transitions. Bottom 
panel: An example H2 transition fitted with the 2 and fiducial 4 component models (dashed and solid 
curves, respectively). Only a small difference is noticeable by eye. Top panel: CRS for the 2, 3 and fiducial 
4 component models (dashed, dotted and solid lines respectively). The significant, many-pixel excursions 
outside the ± 1 -a range for the 2 and 3 component models indicate unmodelled velocity structure, not readily 
apparent when inspecting the transitions "by eye". 



the model with the lowest x 2 P er degree of 
freedom, xl, was selected as the fiducial one. 
The fiducial 4-component model was found to 
have a lower x 2 than the 3-component models, 
which in turn had lower xl values than the 2- 
component model. Fits using 5 and 6 VCs were 
attempted but the additional components were 
deemed statistically unnecessary by vpfit. 

We introduce here the 'composite resid- 
ual spectrum' (CRS) as another, complimen- 
tary tool for discriminating between different 
absorption models. It is constructed by av- 
eraging the 1-cr normalized residual spectra 
(i.e. [data - fit] /error) of regions containing ab- 
sorption lines of interest on a common velocity 
scale. A CRS is a better visual diagnostic of 
poorly fit velocity structure, especially under- 
fitting, than "by eye" inspection of a large num- 
ber of H2 transitions. It can be used to aid the 
construction of new fits. Figure|2]shows a com- 
parison of CRS plots for 2, 3 and 4 compo- 
nent fits. It is clear that the 2 and 3 component 
models fail to reproduce the statistical structure 
of the profile, leaving many-pixel excursions 
outside the expected residual range. Additional 



VCs significantly diminish these deviations, 
with no obvious evidence for unmodelled, sta- 
tistically significant structure in the 4 VC fit. 

The value of A/i//i is simply one of the free 
parameters in the fiducial 4 component absorp- 
tion model for which x 2 has been minimized: 

Afi/iu = (+5.6 ± 5.5 slat ) x 1(T 6 . (2) 

3.3. Simulations ofx 1 minimization 

To verify the uncertainties on A/u/fi returned 
by vpfit we performed a Monte Carlo simu- 
lation of the fiducial 4 component absorption 
model used to derive our constraint. Every ab- 
sorption line and fitting region in our model is 
used as the basis of a simulated spectrum. In 
each of 420 realisations the original flux error 
array was used to generate Gaussian noise in 
the simulated spectrum. Note that the error ar- 
ray was scaled by a factor of 0.8 to increase 
the statistical stability of the fits. The molec- 
ular lines in the simulated spectra use a A///// 
value of 5 x 10~ 6 to verify that vpfit is capable 
of detecting a non-zero deviation in fi from a 
zero starting point. Our fit is one of the largest 
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and most complex constructed in a varying- 
fi study, with a large number of fitted pixels, 
free parameters, various links between param- 
eters, zero and continuum level adjustments 
contributing to the final measurement of A/j//i. 
This is the first time vpfit has been tested in 
this computationally-expensive regime. 

The results are plotted in Figure [3] The 
mean 1-cr uncertainty derived on individual 
Aju//i measurements is 4.4 x 1CT 6 (the expected 
value given a scaling of 0.8 of the flux errors). 
This corresponds well to the standard devia- 
tion of 4.1 x 1(T 6 for the 420 Afi/fi measure- 
ments and is co nsistent with the findings of 
iKing et all (l2009h in that the uncertainties re- 
turned by vpfit may be slightly conservative. 
The mean Afi//j value retrieved is 4.8 x 10" 6 . 



3.4. Potential systematic errors 

A wide range of internal consistency checks 
were run: only relatively 'unblended' transi- 
tions were fit; the contribution of certain sub- 
sets of lines (higher /-levels, HD lines, Werner 
transitions) to Afi/fi was turned off; oscillator 
strengths of all transitions were set to be free 
parameters; separate Afi/fi parameters were fit 
for different /-levels; b parameters were al- 
lowed to very between different /-levels. Some 
of these checks relaxed the assumptions made 
in Section 13.11 We found no evidence to in- 
dicate any obvious problems with our data or 
analysis. Details of these and other tests con - 
ducted will be described in M alec et all d2010h . 

Many H2 and HD transitions are used, over 
a wide wavelength range, so the effect of cal- 
ibration errors in the wavelength scale is sub- 
stantially reduced. However, systematic errors 
due to long- or short-range calibration distor- 
tions are possible and were investigated. 

ThAr lamp exposures were used to cali- 
brate our science ex posures. We u s e the tech- 
nique employed in Mur phy et al.l (l2007h de- 
signed to estimate potential long-range distor- 
tions of the wavelength scale. The possible 
wavelength distortion in the region of molec- 
ular absorption was found to be < 30 ms -1 , 
corresponding to a systematic error in Afi/fi of 
+2.0 x 10~ 6 at most. 
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Fig. 3. 420 Monte Carlo simulations of the fiducial 

4 component absorption model used in our analysis 
- verification that vpfit returns reliable, if not con- 
servative, uncertainties. A variation in \i of A////J = 

5 x 1(T 6 was input into the simulations. The mean 1- 
cr error and the standard deviation (labelled as rms) 
are also shown and are found to be consistent. The 
mean A/i/fi value retrieved (vertical dotted line) cor- 
responds well to the value input in the simulation. 



A distinctly diffe rent type of distorti on was 
recently identified bv lGriest et al.l d2010D in the 
HIRES instrument. Comparisons of ThAr and 
iodine cell exposures revealed that transitions 
at echelle order edges appear at negative ve- 
locities with respect to transitions at the or- 
der centres when calibrated with a ThAr ex- 
posure. A calculation of the possible effect on 
A/j./fj. was obtained by introducing a 'counter- 
distortion', with am plitude and shape similar 
to that described in iGriest et all The value of 
A/j./fi derived from this 'corrected' spectrum 
was (+3.7 + 5.5) x 10" 6 , implying a systematic 
error on Aju/// of approximately +1.9 x 10" 6 . 

Adding these two main systematic error es- 
timates in quadrature, we obtain our final re- 
sult: 



LnlH = (+5.6 + 5.5 stat + 2.7 sys ) x 10" 



4. Conclusions 



(3) 



We have measured A/z/yU from a new H2 ab- 
sorber z a bs = 2.059, with a large number of 
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well-defined H2 transitions and including some 
HD transitions, from a very high resolution 
(R * 110000) Keck/HIRES spectrum. The fi- 
nal value in equation ® includes the formal 
1-cr statistical error and systematic errors es- 
timated from possible long- and short-range 
wavelength calibration errors. This is the first 
Keck-based constraint with precision compa- 
rable to those o btained from VLT spectra by 
iKing et all d2008l) who use a similar technique 
to derive A/u/fi. Our new null constraint is con- 
sistent with the null constrain ts from the VLT 
studie s, in cluding those of [Thompso n et alj 
d2009l) and lWendt & Reimersl d2008h . 

At z < 1, comparison of the radio inver- 
sion transitions of NH3 with molecular rota- 
tional lines has yielded two very strong con- 
straints, Au lu = (0.74+04 7 stat ±0.76 svs )xl0~ 6 
at z = 0.68 (Mu rphy et aill2008l) and (+ 0.08 + 
0.47 sys )xl0" 6 atz = 0.89 dHenkel et alll2009h . 
While these clearly have superior precision 
and smaller potential systematic errors than the 
H2/HD constraints, direct comparison is dif- 
ficult because of the possibility, in principle, 
for spatial fi variations: the different molecu- 
lar species (NH3 and H2/HD) trace regions of 
different densities and, therefore, different spa- 
tial scales and environments. If /v does vary, we 
do not know what that variation depends on, so 
it is presumptuous to prefer one type of mea- 
surement over the other. Clearly, much larger 
samples of both NH3 and H2/HD constraints, 
in overlapping redshift ranges, would allow ad- 
ditional tests for systematic errors and stronger 
conclusions to be drawn. 

It is also important to reduce the statisti- 
cal errors in individual H2/HD measurements 
by substantially increasing the SNR of the 
optical spectra. We have demonstrated here 
that, at least for our spectrum of J2 123 -0050, 
the statistical error dominates over system- 
atic errors. Nevertheless, if such systematic er- 
rors are not also reduced by improving the 
SNR, it becomes more important to increase 
the number of H2/HD absorbers for measur- 
ing Ayu///. Surveys targeting particularly gas- 
and/or metal-rich absorbers have been suc- 
cessfu l in discovering more H2/HD systems 
(e.g. iNoterdaeme et all 120081) but the number 
known in the northern hemisphere is small. 
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